How to find conductance tensors of quantum multi-wire junctions through static 
calculations: application to an interacting Y junction 
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Conductance is related to dynamical correlation functions which can be calculated with time- 
dependent methods. Using boundary conformal field theory, we relate the conductance tensors of 
quantum junctions of multiple wires to static correlation functions in a finite system. We then 
propose a general method for determining the conductance through time-independent calculations 
alone. Applying the method to a Y junction of interacting quantum wires, we numerically verify 
the theoretical prediction for the conductance of the chiral fixed point of the Y junction and then 
calculate the thus far unknown conductance of its M fixed point with the time-independent density 
matrix renormalization group method. 
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Advances in molecular electronics can extend the lim- 
its of device miniaturization to the atomic scales where 
entire electronic circuits are made with molecular build- 
ing blocks [l| . Single molecule junctions connected to two 
macroscopic metallic leads have already been successfully 
fabricated Q . A key ingredient of any such molecular cir- 
cuit is a junction of three or more quantum wires where 
electrical current is conducted through a molecular struc- 
ture between the wires. Determining the conductance of 
such junctions is a long-sought and challenging goal. 

Landauer-Biittiker's formalism does not account for 
electron-electron interactions which play a key role in low 
dimensions. Functional renormalization group methods 
have been helpful in studying the interaction effects Q 
but they are also dependent upon the presence of nonin- 
teracting leads. There are a number of challenging prob- 
lems in junctions of multiple interacting wires. For in- 
stance, a time-reversal symmetric Y junction of Luttingcr 
liquids has a nontrivial fixed point known as the M fixed 
point 0], which has remained an open quantum impu- 
rity problem. For spinfull electrons, there are nontrivial 
fixed points even for junctions of two quantum wires 
Numerical methods like density matrix renormalization 
group (DMRG) 0,0) could potentially be efficient tools 
for computing the conductance of junctions with an ar- 
bitrary number of wires and interactions. Conductance 
is a property of an open quantum system and is related 
to dynamical correlation functions. It may thus appear 
that one needs the more computationally demanding dy- 
namical methods such as time-dependent DMRG for con- 
ductance calculations [i- 14 1. These require large systems 
and averaging over a time window. 

The goal of this Letter is to find generic relations that 
permit the computation of conductance from static equi- 
librium calculations alone, thus making it possible to ob- 
tain the conductance tensors of arbitrary complex junc- 



tions with any number of leads using time-independent 
methods, such as standard DMRG. As a concrete ap- 
plication, we are interested in the open problem of the 
conductance of the M point of a Y junction of quan- 
tum wires. We argue that the T = linear conductance 
Gij (defined by U — JY Gij Vj where Vj is the voltage 
applied to wire j and Ij is the inward directed current 
on it) is determined by the conformally invariant bound- 
ary condition (BC) describing the infrared fixed point of 
quantum junctions. Finding the conductance could then 
be helpful in determining nontrivial fixed points. 

Using boundary conformal field theory, we find a 
generic relationship between the dynamical correlation 
functions of a semi-infinite quantum junction and the 
static correlation functions of a finite system obtained 
by a conformal transformation. This relationship allows 
us to extract certain coefficients from a time-independent 
numerical calculation of the ground state expectation val- 
ues of appropriate local operators which uniquely deter- 
mine the dynamical correlation functions and give us the 
linear-response conductance of the junction through the 
Kubo formula. After establishing the generic continuum 
formalism, we discuss the application of the formalism to 
a discrete lattice computation. 

We present DMRG results on the application of the 
method to an interacting Y junction of quantum wires. 
We successfully verify the theoretical prediction [3] of 



1) 



for the chiral fixed point of a 
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Y junction and numerically calculate the conductance at 
two values of the Luttingcr parameter g for the time- 
reversal symmetric M fixed point. 

Consider M semi-infinite wires connected to a junction 
with arbitrary structure and interactions. We represent 
the system as a quasi-one-dimensional structure with all 
the wires running parallel to the positive x axis (Fig. Q] 
left panel). The wires are connected to the molecular 
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structure of the junction at x — 0. Setting the velocity of 
the charge carriers to unity, we introduce the complex co- 
ordinates z = t + ix, where r is the imaginary time. The 
multiple wires are described by M species of electrons 
living on the upper half complex plane. Notice that spin 
can also be taken into account by doubling the number 
of species, but here we focus on spinless electrons. 

The low energy behavior of the system is given by a 
renormalization group fixed point. Here we treat the 
junction problem as a boundary conformal field the- 
ory [Tif . In the bulk of the wires, we have a conformal 
field theory (CFT) with central charge M and bosonic 
fields <j>j , j = 1, . . . , M obtained from the fcrmionic fields 
via standard bosonization. We assume that the effect 
of the junction (the molecular structure connecting the 
wires) on the fixed-point behavior is generically encoded 
in a conformally invariant boundary condition for the 
CFT in the upper half-plane and represented by a bound- 
ary state \B) [ljjj. This assumption is at the heart of a 
paradigm in quantum impurity problems and has been 
repeatedly verified in a multitude of such problems [I?} • 

Let us first consider the CFT in the absence of the 
boundary. In terms of the bosonic fields <f)j or their dual 
fields 0j, the action for this system is given by 
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where g is the Luttingcr parameter. For nonintcracting 
wires g = 1 and we have g < 1 (g > 1) for repulsive 
(attractive) interactions. The primary operators for this 
CFT are the vertex operators e lo "^ j (the fermion creation 
and annihilation operators are of this form) and the cur- 
rents 

4(z) = -^dP&z) J R (z) = -*8P{z,z) 
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where d = d z — \{d T — id x ) and 8 = d s — \{d T + 
id x ). In the absence of a boundary, the only nonvanishing 
correlation functions are the chiral ones: 
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we have Ag = ^jf%^r^ where \0, 0) is the highest 
weight state corresponding to a generic operator O (here 
O = J^Jfi) Note that at the fixed point the corre- 
lators are determined by conformal symmetry up to the 
coefficients above and as long as we are interested in the 
fixed-point conductance of a given junction, a direct nu- 
merical calculation of dynamical correlation functions is 
not required. The fact that the dynamical correlations 
are directly related to equal-time correlations by confor- 
mal invariance is the key reason why we can obtain the 
conductance from a time-independent calculation. 

If we know the coefficients Ag , we can use the Kubo 
formula [Jl 
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with J 3 ' = J R — J J L to obtain the conductance of the junc- 
tion upon performing the integrations. Let us focus on 
calculating the off-diagonal elements of the conductance 
tensor Gij. The diagonal elements arc not independent 
since Gij = Y) j Gij = due to current conservation 
and the fact that a common voltage applied to all wires 
gives zero current. For i ^ j, as seen in Eq. ([I]), the 
correlators of the same chirality vanish. Using the Kubo 
formula above and after doing the r integral by contour 
integration we obtain for i ^ j 

j-L 
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where H is the Heaviside step function. Since both x and 
y are positive, H(—x — y) = in the second term. 

The next task is to extract A^ by measuring an equal- 
time ground state correlation function. The conductance 
is well-defined for a semi-infinite system, while the equi- 
librium equal-time correlators are easily obtainable in a 
finite system. To get around this issue, we use the fol- 
lowing conformal transformation to map the upper half- 
plane z — t + ix into a strip of width £, 
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and similarly for (J~ T Jft(zi) J 3 R (z2)) where T T indicates 
imaginary time-ordering. 

Let us now consider the system on the upper half- 
plane. In general, the presence of the boundary does not 
change the correlation between primary operators of the 
same chirality [3| but introduces additional correlations 
between the left movers and the right movers, 
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The coefficients A% are determined by the BC on the 
real axis. In fact, in terms of the boundary state |£>), 



Notice that for w = u + iv we have < v < I. The 
mapping takes the positive real axis to the boundary v = 
and the negative real axis to v = £. 

Now consider a primary operator O (here O = J l L J R ) 
with 
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in the semi-infinite plane with BC B on the real axis. 
Using the transformation (0(w)) = \^\~ x ° (0(z)) we 
obtain after some algebra 
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FIG. 1. The conformal mapping from the upper half-plane 
to the strip and the corresponding physical systems. The 
actual junction (here a Y junction, for instance) is illustrated 
next to to the corresponding complex plane regions. The left 
panel shows the semi-infinite junction (the upper half-plane) 
and the right panel shows the finite system obtained by the 
conformal mapping (the strip). 



Since the zero temperature finite system is invariant un- 
der translations in it, the static ground state expectation 
value of the local operator J l L {x)J J R (x) with x the dis- 
tance from a boundary in the finite system of length £ is 
expected to behave as --jfy-A^ [2 sin ( fx) / f] '. This 
ground state expectation value permits the determina- 
tion of the A l g coefficients and thus the conductance via 
Eq. ([2]) using solely static equilibrium computations. 

In the remainder of this Letter, we discuss the appli- 
cation of the method to a tight-binding lattice numerical 
calculation. First we consider the BCs for the finite sys- 
tem obtained from the transformation Eq. ([3]). Because 
we are assuming a conformally invariant BC, the bound- 
ary v = will have the same BC B as the real axis of the 
semi-infinite plane. To determine the BC for the v = £ 
boundary created by the transformation, we consider the 
fcrmion creation and annihilation operators which trans- 
form as 
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Using m = 42 n = ie" 
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_i. e -f« j we fj nc j that up to a sign coming from choos- 
ing the branch of the square root in the transformation 
above, the BC at v = £ is the same as v — 0. This 
does not mean however that in a microscopic implemen- 
tation, one should place the mirror image of the junction 
at v = i. Notice that the chiral fields switch role for the 
two boundaries; i.e. the right-movers are the outgoing 
(incoming) states for the v = (v = £) boundary. So for 
instance if the BC at v = is described by an S-matrix 
S, the BC at v = I is given by the inverse S-matrix S' 
up to a sign. 

We consider a half-filled system with an even number 
of sites in each wire and particle-hole symmetry in the 
bulk of the wires. The Hamiltonian of the system is given 
by Hl + Hb + Hr where Hr is the bulk Hamiltonian of 
the wires and Hl.r are the junction Hamiltonians. In 
terms of the fermionic operators ^l 3 on lead j and site 



where n J m = 'J'- 7 "F^ , the hopping amplitude is set to 
unity and V is the interaction strength. To implement 
the correct boundary condition at v = £, we construct Hr 
by applying a particle-hole and a time-reversal transfor- 
mation to Hl ■ As a concrete example we have the bound- 
ary contribution H R = -^ffVij^j \iH L = \\>\ I , for 
a Hcrmitian matrix T. 

The coefficients A 1 ^ can be extracted from the ground 
state correlators of chiral operators in a finite system. We 
cannot directly model chiral fermionic creation and an- 
nihilation operators, but indeed we can model chiral cur- 
rents using density and current operators. The (nonchi- 
ral) current operator on the link between sites m and 
m+ 1 is given by P m = i{^l n+1 ^ m - ^L^ n+1 ). Also, 
for a bond between sites m and m + 1 , the charge den- 
sity operator (with the background charge subtracted) is 

Chi- 
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ral current operators on the lattice are related to the 
density and nonchiral current through 

J*(x) = v(J R (x)-4(x)), N*(x) = J R (x)+4(x) (4) 

where v is the velocity of charge carriers, and x takes 
values in the middle of the bonds, x = m + 1/2. 

The Luttinger parameter g and the velocity v depend 
on the strength of the interaction [2l| . We have from 
the Bethc ansatz 
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arccos (V/2) ' V ; 



We would like to find the off-diagonal elements of the 
conductance tensor Gy. The chiral correlators are pro- 
portional to 5ij in the half-plane and therefore vanish for 
i 7^ j. They also vanish in the finite system since it is 
obtained from the half-plane by a conformal mapping. 
Using Eq. (j4]) , we can then write for i 7^ j 



(J\x).P{x)) = -v 2 {(JUx)J R (x)) + (J R (x)Jl(x)) 
{N*(x)ji(x)) = v ((Ji(x)J R (x)) - (J R (x)4(x))) . (6) 



By measuring the above two correlation functions (the 
second correlator (N l (x) J 3 (x)) is identically zero in the 
absence of a magnetic flux) , we can find the chiral corre- 
lators we need for calculating the conductance. We fit the 
measured {J l L (x)J R {x)) with A [2 sin (fx)/j] to find 
the coefficient A and then obtain the conductance using 
Eq. [21 Note that the fixed point behavior cannot be ob- 
served close to the boundary so the numerical results for 
small x must be excluded in fitting the data. 
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Below we consider an interacting Y junction of three 
quantum wires. As a check, we first apply the method to 
the chiral fixed point whose conductance has been pre- 
dicted in Ref . Q ■ We then present the first numerical cal- 
culation of the conductance of the time-reversal symmet- 
ric M fixed point whose properties have so far remained 
unknown. As seen in the inset of Fig. [2j the junction 
simply connects the three endpoints of the wires with a 
hopping amplitude t (the amplitude is 1 in the bulk) . As- 
suming that the loop formed at the junction is threaded 
by a magnetic flux (f>, the system will renormalize to the 
chiral fixed point for nonzero <f> and 1 < g < 3 indepen- 
dently of the hopping amplitude t. In the same range of 
<7, the system rcnormalizes to the M fixed point in the 
absence of a magnetic flux (0 = 0). 

The numerical calculations for both fixed points were 
performed for g = 1.5, 2.0 and t = 0.7, 1.0. For the chiral 
fixed point, using tfi = ir/2, £ = 39 (40 sites on each 
chain) and the number of DMRG states m = 600, we 
found strong agreement between the numerics and the 
predicted conductance G12 = —2 g^p-(.9 + 1)77 Q as 
seen in Fig. [2] The correlation functions are independent 
of the hopping amplitude t indicating convergence to the 
universal conductance of the fixed point. For the M fixed 
point, using (f> = 0, £ = 59 and m = 1000, we obtain a 
conductance of G 12 = -0.5293^ (G 12 = -0.6146^) for 
g = 1.5 (g = 2.0). None of the properties of the M 
fixed point were previously known in the literature fl7j ]. 
These values of the conductance depend on g nontrivially 
and are different from those of all the other fixed points 
which have been identified so far. Thus our method and 
the numerical results provide evidence that the M fixed 
point is of a different nature. 

In summary, we presented in this Letter a generic rela- 
tionship between the elements of the conductance tensor 
of quantum junctions at a renormalization group fixed 
point and coefficients in the static correlation functions 
of chiral current operators. The method hinges on a fun- 
damental relationship between static and dynamic cor- 
relation functions in boundary conformal field theory. 
We proposed a method for obtaining conductances from 
static calculations by extracting these coefficients using 
nonchiral current and density correlations in a lattice 
computation. We applied the method to an interacting 
Y junction of Luttingcr liquids using time-independent 
DMRG and presented a numerical verification of the the- 
oretically predicted conductance of the chiral fixed point 
as well as a calculation of the previously unknown con- 
ductance of the M fixed point. We emphasize that this 
general method turns the time-independent DMRG into 
a powerful tool for extracting conductances of completely 
generic molecular electronic junctions which constitute a 
wide class of quantum impurity problems. 
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FIG. 2. The inset shows the lattice Y junction. The thick 
bonds at the junction have a hopping amplitude t. The lines 
drawn through the data points show the CFT prediction with 
the exact analytical conductance in the case of the chiral fixed 
point (<f) 7^ 0) and the fit used for obtaining the conductance 
in the case of the M fixed point ((f) — 0). As expected, the 
conductance is independent of t. 
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